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ABSTRACT 

Merging binaries of compact relativistic objects (neutron stars and black holes) 
are thought to be progenitors of short gamma-ray bursts and sources of gravitational 
waves, hence their study is of great importance for astrophysics. Because of the strong 
magnetic field of one or both binary members and high orbital frequencies, these binaries 
are strong sources of energy in the form of Poynting flux (e.g., magnetic-field-dominated 
outfiows, relativistic leptonic winds, electromagnetic and plasma waves). The steady 
injection of energy by the binary forms a bubble (or a cavity) filled with matter with the 
relativistic equation of state, which pushes on the surrounding plasma and can drive 
a shock wave in it. Unlike the Sedov-von Neumann- Taylor blast wave solution for a 
point-like explosion, the shock wave here is continuously driven by the ever-increasing 
pressure inside the bubble. We calculate from the first principles the dynamics and 
evolution of the bubble and the shock surrounding it and predict that such systems can 
be observed as radio sources a few hours before and after the merger. At much later 
times, the shock is expected to settle onto the Sedov-von Neumann- Taylor solution, 
thus resembling an explosion. 

Subject headings: ISM: bubbles; stars: neutron; binaries: close; shock waves; ISM: jets 
and outflows 



1. Introduction 

Neutron star-neutron star (NS-NS) and neutron star-black hole (NS-BH) binaries are of great 
interest in astrophysics. First, they are thought to be progenitors of short gamma-ray bursts 
(GRBs), though the evidence is based mostly on the shortness of the time-scale of the final merger 



process, in addition to the energetics considerations and population studies (Blinnikov, et al. 1984 
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Eichler, et al. 1989, Nakar 2007 Berger 2011). Hence, the natural question arises: Are there any 



observational signatures that can help to tell apart the merging binary progenitor from alternative 
GRB models? Unlike other short GRB studies (including those studying their radio emission, e.g.. 



Pshirkov Sz Postnov 2010; Nakar & Piran 2011), we investigate here what happens to the system 



long before the merger, not during or after. Second, merging binaries are also considered to be 
strong sources of gravitational waves, which can potentially be observed with upcoming detectors 
(e.g.. Advanced LIGO). Cross-correlating the gravitational wave signal with astrophysical objects 
observed using other techniques will advance astrophysics drastically. 

The binary evolution and its orbital period are determined by the emission of gravitational 
waves. Besides, neutron stars possess large magnetic dipole moments, hence the electromagnetic 
energy is also extracted from the system via Poynting flux. Theoretical estimates and numerical 
modeling of NS mergers suggest significant amounts of energy, ~ 10^^ erg ( McWilliams & Levin 



2011 Piro 2012 Etienne, et al. 2012), to be released in the electromagnetic form, and even more 



for magnetar binaries. The actual energy extraction mechanism may be rather complicated. Here 
we make an analogy with pulsar electrodynamics. A rapidly spinning neutron star magnetosphere 
produces a magnetized, relativistic electron-positron wind that exerts spin-down torque and ex- 
tracts rotational energy dGoldreich fc Juhan|p69| ). The presence of relativistic leptonic plasma in 
the magnetosphere changes its structure drastically, which complicates analytical analysis. How- 
ever, despite the structure of the force-free magnetosphere is different from the dipole field, the 
electromagnetic energy extraction rate given by the expression for the magnetic dipole radiation 



in vacuum (Gunn & Ostriker 1969) is only within a factor of two, depending on geometry, of the 



exact rate obtained from direct numerical simulations (Spitkovsky 2006). Thus, for the case of a 



NS binary, we also expect that the electromagnetic energy is extracted in the form of the Poynting 
flux and it is safe to assume that the energy extraction rate — hereafter referred to as the 'Poynting 
luminosity' or just the 'luminosity', L[t) = dE/dt — is of the order of the electromagnetic losses 
in vacuum. Hence a bubble (or a cavity) forms in the surrounding medium. For simplicity, we 
assume spherical symmetry of the bubble and the uniform density of the ambient medium. These 
assumptions can readily be generalized to include environmental effects; for example, rapidly spin- 
ning NSs can themselves produce outflows and form pulsar wind nebulae around them. However, 
the uniform density assumption is reasonable if the neutron stars are non-rotating and do not have 
pulsar wind nebulae around them. The matter composition (relativistic plasma, magnetic field, 
waves) in the bubble depends on the exact mechanism of the energy extraction (e.g., a magnetized 
outflow, relativistic wind, electromagnetic and/or plasma wave emission and/or conversions, etc.), 
as well as the structure, evolution and interaction of magnetospheres of orbiting (and spinning) 
NSs — the problem yet to be solved. However, regardless of the composition, the material in 
the bubble has a relativistic equation of state, 7 = 4/3. Because of the electromagnetic nature 
of the process producing the bubble and because the equation of state of matter inside it is that 
of electromagnetic radiation, we colloquially refer to them as the 'electromagnetic bubbles'. Being 
a strong function of the orbital frequency, the binary Poynting luminosity L(t) and, hence, the 
pressure inside the bubble are rapidly increasing around the merger time. The expanding bubble 
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Fig. 1. — Schematic representation of the bubble+shock system around a NS-NS binary. 

pushes onto the ambient medium at an increasing rate and can drive a shock in it, as is ihustrated 
in Fig. [T] As the shock forms, it can be detected, e.g., via synchrotron radiation from accelerated 
electrons in shock-amplified magnetic fields. 

In this paper we calculate the dynamics of the electromagnetic bubble and the bubble-driven 
shock around merging double neutron star or magnetar binaries and make observational predictions. 
In Section [2] we derive the evolution of the binary orbital period assuming the orbits are circular. 
In Section [3] we discuss various processes of electromagnetic energy extraction and evaluate the 
Poynting luminosity. Section |4] presents the model of the bubble+shock model and gives its evo- 
lution in the analytical form, including the finite-time-singular solution. Sections [5] and |6] provide 
numerical estimates of the bubble+shock system parameters and make observational predictions, 
respectively. Section [7] presents discussion and conclusions. 



2. Binary inspiral 

The dynamics and merger of a NS binary is determined by emission of gravitational waves. 
Here we assume circular orbits of the NSs, for simplicity. We also neglect general relativistic effects, 
though they will be important in the final few seconds of the merger. We also consider the NSs 
as point masses, hence neglecting tidal effects. With these assumptions, the energy loss rate in 
the systems of two gravitating bodies of masses Mi and M2 orbiting in circular orbits about their 
common center of mass is 

— = = ^m^r^Vt^. (1) 
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Here m = MiM2/(Mi + M2) is the reduced mass, is the orbital frequency, G is the gravitational 
constant, Qao- = fn{3xxXfj — r'^Sx(^) is the quadruple moment, r = ri — r2 is the mass separation vec- 

* ■ * 2 * ■ * X(T * ■ ** 

tor with r and x\ being its magnitude and components, Q = Q Q^^ and (Jao- is the Kronecker ten- 
sor. By the virial theorem mv'^/2 = K = -E = -U/2 = GMiM2/2r = (m/2) [G(Mi + Ma)^]^/^. 
NSs are relativistic objects so it is convenient to normalize the distance and the frequency by the 
radius of the maximally spinning Kerr black hole having the total mass of the system, Rg, and the 
Keplerian frequency at Rg-. 

Rg = G{Mi + M2)/c\ Qg = [G{Mi + M2)/RI] . (2) 

Then 

r/Rg = (n/ng)-^/^ and w/c = (3) 

In the i?g-units, the radius of a neutron star is Rns ~ '^Rg (for equal mass stars) and the shortest 
separation distance in the binary (when NSs merge) is Rm ~ 2Rns ~ QRg- Realistically, disintegra- 
tion of the neutron stars may occur at even larger separations due to strong tidal forces. Moreover, 
general relativistic effects become important at the final stage of the inspiral and merger, but they 
are omitted from consideration in our simple model. We parameterize the shortest separation 
distance with k/j as 

Rm = K-nRg, KR ~ 6. (4) 

Note that the smallest separation between the solar mass NS and BH in the NS-BH binary is 
smaller, Rm ~ Rns + Rbh ~ i.e., kr ~ 4. 



From equation noting that Rg^g = c and —E = K = {mc? /2){VL/Vtg)'^/^ , one has the 
equation for the orbital frequency 



d{9./Q.g) 96 M1M2 ^ ( ^ 



11/3 



7^ ■ (5) 



dt 5 {M1 + M2Y " ^ V^s 

It has a general solution with the finite time singularity 

VL{t) = ^i{l-t/ts)-'"\ (6) 

where Q,i is the initial orbital frequency and the source time is 

5 {M1 + M2Y fn,^'^^^ 



(7) 

256 M1M2 yngj 3 ■ ^ ^ 

GRg there is maximum frequency and time beyond which equation 
(|6]) is inapplicable, 

Q,m/^g = ^ , tm/ts = 1 — I'^'fl i^i/^g)^^'^ ■ (8) 
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3. Poynting luminosity L{t) 

A merging NS binary is a system of orbiting (and possibly spinning) magnetic dipoles. Be- 
cause of strong magnetic fields in and rapid motion of the magnetospheres, the induced electric 
fields are strong enough to develop an electron-position cascade ( Timokhiri||2010 ) , which loads the 



magnetospheres with leptonic relativistic plasma; hence the force-free configurations are formed and 
maintained. Therefore, the electromagnetic (EM) energy extraction can, in general, be different 
from the simple dipole/quadrupole radiation losses. Force- free modeling of magnetospheres of close 
NS-NS has not been done, so the details of the EM energy extraction remain unclear. However, we 
can use the analogy with pulsar electrodynamics, which indicates that the energy is extracted in the 
form of relativistic leptonic winds and field-aligned currents — generally referred to as the Poynting 
fiux — and, moreover, the EM losses are only a factor of two or less away from the low-frequency 



electromagnetic wave emission losses due to the magnetic dipole radiation in vacuum (Spitkovsky 



2006), though some studies (Timokhin 2006) suggest a more complicated picture. Therefore, we 
use the standard EM losses (dipole, quadrupole) in our estimates of the power taken away via the 
Poynting fiux. One should keep in mind that these are estimates only; the actual mechanism and 
the electromagnetic energy extraction rate are unknown at present. 

In general, the emitted power depends on the orientation of the dipole moment, the orbit 
eccentricities and the angular momentum of the system, as well as the spins of the neutron stars. 
Below, we consider the induced electric dipole and magnetic dipole radiation losses, and we briefiy 
discuss quadrupole radiation losses. If the NSs are not rapidly spinning, the most efficient EM 
energy loss is the induced electric dipole emission, which occurs due to the orbital motion of 
magnetic moments of NSs. It is present even if the NSs are non-spinning at all. However, the 
electric dipole is of the order of ~ v/c of the magnetic moment. Hence, the electric dipole emission 
contribution can be sub-dominant if the NSs are rapidly spinning. Moreover, it is natural to expect 
that during the final moments of binary inspiral, the NSs may be tidally locked and their angular 
frequency will be equal to the orbital angular velocity. In this regime, the magnetic dipole emission 
losses will dominate and the induced dipole emission will be jc? times weaker. Quadrupole 
emission losses will also be present, but it is naturally weaker by a factor of jc? and effectively 
renormalizes the induced electric dipole emission power. 

It will be shown that in all cases, the Poynting luminosity of the binary obeys the finite-time- 
singular (FTS) law, cf., Eq. 

L(t) = L,(l-t/t,)-P, (9) 

where the index -p depends on the type of radiation mechanism, Lg is the initial luminosity (at 
t = 0) and ts is the binary lifetime - the time until merger in the Newtonian approximation with 
the binary members being point-like objects. In reality, they are of finite size, so the actual merger 
time tm is less than tg, as we discussed in the previous section. From equations ([T]) and ([s]), one 

256 M1M2 ^ ' 
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Note that this quantity is independent of initial conditions of the binary insipral, i.e., of the initial 
orbital frequency, Qi. The total energy output does not diverge and is estimated to be 

Es = / Ls{l-t/ts)-Pdt 
Jo 

(p-l)-iLA(l-Wis)-(^-'\ (11) 
where the last approximate equality holds true if p > 1 and ts — tm tg- 



3.1. Induced electric dipole losses 

Here we consider the simplest case of non-spinning NSs with magnetic moments and fj,2 
aligned with the orbital angular momentum. A magnetic moment moving non-relativistically with 
a velocity v induces an electric dipole moment d = (v^^ x /x)/c. For the NS system at hand, 
di = rifiivi/c and d2 = r2/i2?^2/c, where vi = ri = rM2/(Mi + M2), V2 = r2 = rMi/(Mi + M2) 
and "hat" denotes a unit vector. 

The dipole radiation power is 

where we used that d = {fiv/c)e^^^. The Poynting luminosity of the NS binary is the sum of the 
luminosities of each of the dipoles. 



_ 2 (MfMl + ^jMf) ^ ^14/3 

'^^'>- 3c^ (Mi + M2)2 [nj ' 



(13) 



where we used Eq. Q. With equation ([6]), the luminosity becomes 

L(^) = L,(l-^/^.)-'/^ (14) 

where 

_ 2 {JiMl + ^4Mll ,fn,Y/' 

^^-3c3 (Mi+M2)2 ^^^UJ ■ ^ ^ 

Using Eq. ([s]) and that Rg^g = c, one can estimate the peak luminosity = L{tm) to be 



ng 



2 ng {niMi + ^^Mi^ 
3^^ {M1 + M2Y 



(16) 
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Note that this quantity is independent of the initial orbital frequency, j . The total energy output 
is estimated to be 

192 ng M1M2 ' ^ ' 

We note that in this case the Poynting luminosity index is p = 7/4 in Eq. 



3.2. Magnetic dipole losses 

We can also assume that the neutron stars are tidally locked and their magnetospheres rotate 
with the orbital frequency and result in the magnetic dipole losses with the power 

dE^M 2 ..2 2 /io\ 

Using equation ([6]), one obtains the luminosity of two tidally locked NSs 

L(t) = L,(l-t/t,)-3/2, (19) 

where 

= j3(/^? + ^i)^^(|)'- (20) 

The peak luminosity = L{tm) is obtained using Eq. Q to be 

- ' ^H^^^+f^^). (21) 



Finally, the total energy output is 

^ 54 (Mi + M2)^ 

128 M1M2 ■ ^ ' 

and we remind that the magnetic dipole losses result in the Poynting luminosity index is p = 3/2 
in Eq. 



3.3. Magnetic quadrupole losses 

An orbiting magnetic dipole has a magnetic quadrupole moment Q ^ fir. Quadrupole losses 
are generally weak, but the induced electric dipole losses are, in fact, of the same order mv/c. Here 
we just look for a scaling. The emitted power is 

^= 1 QVr2l76ocf]i4/3, (23) 
at 45c'' 

which yields p = 7/4 as in the induced electric dipole case. Thus, this energy loss channel just 
makes a correction to the induced electric dipole one. 
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3.4. Comments on magnetar binaries, NS-BH binaries some NS-NS binaries 

Magnetars are neutron stars with the surface magnetic fields of the order of or larger than 
the Schwinger field Bs — 4:A x 10^'^ G. All the estimates in previous sections remain true with 
the magnetic moment being about three orders of magnitude larger than for normal NSs, ^Mag ~ 
1033 G cm . This increases the luminosity and the total energetics by a factor of one million, 
E (X fi^ and hence E ~ 10^'' erg, as is estimated in Section [sj 

Evolution of a NS-BH binary is very similar to that of a NS-NS binary, except for the last 
moments before their merger when GR effects are important, hence equation ([g]) holds. Poynting 
luminosity produced by a NS-BH binary is mostly via the electric dipole losses due to the electric 



dipole induced on a BH horizon, so one has (see Eq. (5) of McWilliams Sz Levin 2011) 



L oc v^r-^ oc n^^/^, (24) 

which also yields p = 7/4. However, if the NS companion is tidally locked, then the magnetic dipole 
losses would dominate the Poynting luminosity. 

Finally, in NS-NS binaries, one of the members can have a substantially weaker magnetic field 
then another, as is the case in the double pulsar system PSR J0737-3039, for example. In this case, 
the stronger- field binary member induces the electric dipole moment on the other NS, very much 



like in the NS-BH binary. For such a case, a unipolar inductor model has been proposed (Piro 



2012). This model generally predicts the Poynting flux luminosity to be similar to the case of the 



electric dipole losses, Eq (13). 



4. Bubble-|-shock system evolution 

The formation and dynamics of bubbles and bubble-driven shocks has been studied in detail 

. Here we briefly outline the idea and use appropriate 
results. For simplicity, we assume that the electromagnetic bubble is spherical and expands into the 
interstellar medium (ISM) of constant density, as shown in Fig. [l| which is a good approximation 
for non-spinning neutron stars which do not have pulsar wind nebulae. The bubble is fllled with 
matter with the relativistic equation of state parameterized by the adiabatic index 7em = 4/3. The 
bubble surface acts as a piston and exerts pressure on the external medium producing an outgoing 
strong shock [the compression ratio is k = (7 + l)/(7 — 1) ~ 4] in the cold unmagnetized ISM with 
mass density pisM and the adiabatic index 7 = 5/3. The shell of shocked ISM is located in between 
the shock and the bubble. The mass density and pressure of the gas in the shell are determined 
by the Rankine-Hugoniot jump conditions. The pressure equilibrium throughout the system and 
across the bubble-shell interface (i.e., a contact discontinuity) is assumed. 

The time-dependent Poynting power L{t) goes into the following components: the internal 
energies of the bubble, dC/bubbie, and the shocked gas shell, dUsheii, the change of the kinetic energy 



in a separate paper (Medvedev Sz Loeb||2012 
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of the bulk motion of the shell, dKshcW, assuming its swept-up mass Mg^ept is constant, as well as 
heating, dU@shock, and acceleration, dK^shock, of the newly swept ISM gas at the shock. The pdV 
work due to the expansion can be neglected because the external pressure is vanishing in the cold 
ISM. Thus, the master equation is 



L{t) dt = dC/bubbie + dUsheW + dK, 



shell I M= 



wept 



+ dU<r 



'shock 



+ dK, 



Ishock' 



(25) 



All the quantities can be expresses as a function of one dependent variable, the shock radius R{t), 
for example. Other quantities follow straightforwardly, e.g., the shock velocity is v = R, the bubble 



radius Kb = {1 



physically transparent: L{t) 



(3/4) ^/'^i?, and so on. The structure of the master equation is 
k ~ dtipR^v'^) ~ pisuiciR'^R^ + C2R^RR), where ci and C2 are 



some constants. The actual calculation (Medvedev &: Loeb 2012) yields: 

L{t) 



ciR^R^ + C2R^RR, 



(47r/3)pisM 



where 



6 



ci 



C2 



7em + 1 



+ 2 



(7 + 1)2 \7EM - 1 

4(7EM + 1) 12 
(7 + l)2(7EM-l) 7^-1 



7.6, 



7 + 1 



1/3 



1 



~ 4.6. 



(26) 

(27) 
(28) 



Solution of this inhomogeneous second-order nonlinear differential equation yields the shock radius 
as a function of time for any given luminosity law L{t). 

In Section [3] we have shown that the Poynting luminosity of a NS binary is represented by the 
FTS law, Eq. 

L{t) = Ls{M/t,)~P, (29) 

where Ai = — t. Approximate analytic solutions of Eq. (26) exist at both early and late times 
( Medvedev &: Loebl|2012 ). At early times, t -^tg, the binary has approximately constant Poynting 
luminosity, so the shock radius and velocity are: 



R{t) 
vit) 



Rs {t/tsf^ , 



where 



R. 



3 L.t^V^' 



Rs/ts 



(30) 
(31) 

(32) 



,47r pismJ 

and some numerical factors of order unity are suppressed for clarity. At late times, i.e., around the 
merger time t ~ is, the luminosity increases rapidly, so does the shock velocity, whereas the shock 
radius approaches a constant: 



R{At) 
v{At) 



R,{At/t 
VsiAt/t 





Sj ! 



-{p-l)/2 



(33) 
(34) 
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Note that the early-time solution is self-similar and the late-time one has a finite time singularity 
and therefore can break down if the velocity approaches the speed of light. Formally, it also breaks 
down at small At < i?/c because of the finite time needed for pressure equilibration throughout 
the system, where c is the speed of light (recall, the material inside the bubble has a relativistic 
equation of state). Nevertheless, it can still be approximately true if tg is corrected for the finite 
light travel time, ts ^ ts + R/c. 



5. Numerical estimates 

In previous sections we made theoretical estimates of the binary evolution, its Poynting lumi- 
nosity and the evolution of the bubble+shock system. Here we make order of magnitude estimates. 

First of all, if the masses of the compact companions are Mi = M2 = tuMq, then, from Eq. 
([2]) one has 

i?g ~ 3 X 10^ m cm, ~ 10^ m'^ rad s"\ (35) 
The binary lifetime, Eq. ([T]), depends on the initial orbital angular speed, Qi, 

t, ~ 2 X 10^ m-^/^j^r^/^ s (36) 

which is defined by the moment when — ?■ 00 or the separation of two point-like masses approaches 



zero. In nature, because of the finite size of the objects, the actual merger time tm, Eq. ( 10 ), occurs 
before tg'. 

Atm = ts-tm^ 10"^ m s. (37) 
The maximum orbital frequency, Eq. ([s]), is 

0^ ~ 7 X 10^ rad s"^ (38) 

The orbital frequency evolution, Eq. ([6]), can be re- written as 

n{At) = nm{At/Atm)-^/^ (39) 
~ 5 X 10^m"^/^At"3/^ rad s~\ (40) 

where At = tg — t and other quantities are in CGS units unless stated otherwise. 

Next, we estimate the Poynting luminosities and bubble+shock parameters for two cases, 
depending on the dominant mechanism of the electromagnetic energy extraction: the magnetic 
dipole and induced electric dipole losses. These two cases correspond to different physical scenarios. 
Case 1 represents the binary in which the NS spin periods are approximately equal to the orbital 



period, Qns ~ ^- Bildsten Sz Cutler (1992) demonstrated that tidal locking is impossible in NS-BH 
binaries and is unlikely (though cannot be completely ruled out due to our ignorance about NS 
internal viscosity) in NS-NS binaries. However, complete locking, H-ns = is not required in this 
scenario. Partial synchronization will lower Poynting losses somewhat and can alter the temporal 
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index, only. Rather strong partial synchronization, (Jl — O^vs)/^ ~ 10% has been estimated in (Piro 
2012) for the binary with orbital frequency J7 ~ 10^ s~^ in which one of the NSs is non-magnetized 



Hence Case 1 is entirely viable, and it is even more so given the lack of detailed understanding of 
interaction of NS magnetospheres. Case 2 corresponds to several physical scenarios: (i) the binary 
with slowly-spinning or non-spinning neutron stars, (ii) the NS-BH binary and (iii) the NS-NS 
binary with one of the NSs having much weaker magnetic field then the other. We note that Case 
1 is more energetically efficient and effectively represents the upper bound on the process, whereas 
Case 2 is presumably more realistic. 



5.1. Case 1: magnetic dipole radiation losses 

In this scenario, neutron stars are approximately synchronized, i^NS ~ hence the magnetic 
dipole losses dominate Poynting luminosity. Since this mechanism is the most energetically effective. 
Case 1 represents the order-of-magnitude upper limit on electromagnetic processes in the binary. 



The electromagnetic luminosity, Eq. (19), is 

L{t) = Ls{l - t/ts)'^/^ = L,n (At/At„)-3/2 _ (4^) 



From Eq. (21), the maximum luminosity, which occurs at the time of merger L{tm), is 

Lm,NS ~ lO^^/xiom"^ erg s"^ (42) 

for a NS-NS binary and an order of magnitude larger for a NS-BH binary (because of lower kr). 
Here we assumed that the NS magnetic moments are = /i2 = M ^^^d have a typical value of 
= fJ'SoW^^ G cm^. For a NS-magnetar or a double-magnetar binary, the luminosity is orders of 
magnitude larger: 

LmMn^g ~ lO^^/^is™"^ erg s"^ (43) 
for a nominal value of the magnetic moment of /^Mag ~ 10^^ G cm^. 



The total EM energy produced in the process, Eq. (22), is 



Es,NS ~ 2 X lO^Vao"^"^ erg, (44) 

i.e., Es^NS ~ 10^^ erg for a typical NS binary, but can be as large as -Bs^Mag ~ 10^'' erg for a 
magnetar binary. We stress that all calculations are done in the non-relativistic limit and do not 
account for general relativistic effects. 



The normalization Lg in Eq. (41) depends on the binary lifetime, which is uncertain in reality. 
In fact, it is unreasonable to take the time since the binary was formed, because the Poynting 
luminosity is very low and other processes determine its ambient medium conditions. For example, 
motion of a binary in the ISM may disrupt and destroy the bubble by ram pressure if the bubble 
expansion velocity (at any moment of its evolution) is smaller than the bulk motion of the binary 



- 12 - 



as a whole, which can be assumed to be a few tens km/s. As we will see below, this condition is 
satisfied if ts is about tens of years. For concreteness, we assume here tg ~ 10 years. One obtains 



)3/2 
-5/2^-3/2 



6 X 10^6 ^iom--"t,i5y ergs 



(45) 



where ts,ioy = ts/{W yr). 



The shock radius and velocity evolve according to Eqs. (30), (31) at early times t <C and 
according to Eqs. (33), (34) at later times around the merger time, t ^ ts- The characteristic 
values are given by Eqs. (32): 



Rs 

Vs 



1/2 -1/5 ,3/10 



inlS 2/5 

iU /igQ m "'ISM,0''s,10j/ 

4 X 10*^ f4o^rn- 



1/2 -1/5 ,-7/10 -1 
"'ISM,0'-s,10y ^™ S ' 



(46) 
(47) 



that is, a typical size of the shock is Rs ~ 70 AU and its minimum velocity Vs ~ 40 km/s for the 
assumed parameters. At late times, the shock scalings with time are 



R{At) = Rs, v{At) = Vs {At/ts 



-1/4 



(48) 



5.2. Case 2: Induced electric dipole losses 

In this case the induced electric dipole (together with magnetic quadrupole) mechanism dom- 
inates, hence 

L{t) = Ls{l - t/ts)-"^ = L„ {At/At^)-"^ . (49) 



The maximum luminosity, Eq. (16), is 



Lm,NS ~ 5 X 10^V30"i ^ erg s ^ (50) 



for a NS-NS binary and about 5 x lO'*® erg s for magnetars. The total energy, Eq. (17), is 

Es^NS ~ 6 X lO^V^o"^"^ erg, (51) 
i.e., Es^NS ~ 10^^ erg for a typical NS binary, but can reach ~ 10^^ erg for the magnetar case. 

The time ts should be smaller in this case because of the lower overall energetics (see discussion 
in previous subsection). We use ~ 3 years. The normalization of the luminosity Lg in Eq. (49) 



IS 



T — T f-'^/^(f -f "17/4 

~ 2 X 10^^ f4om-"H:l/; erg s'^ (52) 
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The characteristic values of the shock, Eqs. (32), are 



P ^ inl4, 2/5^-9/20^-1/5 ,1/4 

Q in6 2/5 -1/2 -1/5,-7/10 -l 

Vs ~ 3 X 10 ^3(, m / njgj^ gt, !^ cm s , 



(53) 
(54) 



i.e., the bubble and shock are an order of magnitude smaller in this scenario. At late times, the 
shock parameters scale as 



R{At) = Rs, v{At) = Vs {At/ts 



-3/8 



(55) 



6. Observational predictions 



Shocks can be observed via synchrotron radiation produced by shock-accelerated electrons in 
generated or amplified magnetic fields. We assume that the electrons and magnetic fields carry 
fractions eg and of the internal energy density, u^heW ~ PiSMv"^, of the shocked ga^ 



(56) 



where ^eJemeC^ is the average energy of an electron, /3e = Ve/c is the average dimensionless electron 
speed, and ng^sheii = n ?^e,iSM — k nisM is the number density of electrons in the shocked gas shell, 
and K ~ 4 is the shock compression ratio. 

Here we consider only the first scenario with the magnetic dipole Poynting luminosity, since it 
provides the most interesting observations limits. From the observational point of view, the value of 
ts is nearly impossible to determine, whereas the bubble or shock radius can be measurable either 
directly (if the image is resolved) or indirectly (by time variability, for instance). Eliminating ts 



between Eqs. (46), (47) and using Eqs. (48), ([37j) we can express the shock speed via the shock 
size: 

v^7xW^ fi3im-'>/\-^i\R;%^At~^^^ cm s-\ (57) 



where we evaluated the shock speed 10^ seconds before the merger and the neutron star surface 
field is ~ 10^^ gauss. In estimates below we will use this expression with explicit dependence on 



the shock size. If desirable, the dependence on ts can readily be restored using Eq. (46). 



The average Lorentz factor of accelerated electrons is obtained from 

l3ele{At) ~ eeK{mp/me) (v/cf 

~ 4 ee/xilm-5/2nf4_oi^-?5A^-'/^ 



(58) 



^The exact calculation of Wsheii (Eq. 30 in 
suppress hereafter. 



Medvedev & Loeb 



20121, differs by a factor 2/(7^ 



9/8 which we 
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so the bulk of the electrons is mildly relativistic for a nominal ~ 0.5. The sub-equipartition 
magnetic field strength is 

B{At) ~ {SireBmpnisu)^^'^ V 

~ 5 X 10-3 e]i'fi^im-'/%s'M,oK%'W^' G, (59) 

that is B is of the order of 0.1 milligauss for a nominal ~ 10^^ which means the field must be 
generated or amplified at the shock by an instability (such as Weibel, Bell, firehose, cyclotron), 
preexisting MHD turbulence or via other mechanism. 

Relativistic electrons emit synchrotron radiation with the peak of the spectrum being at the 
frequency 

Us{At) ~ {27T)-^f^{eB/mec) 

~ 2 X 10^ elef4irn-''^'n^^,o<il^'^tf/' Hz, (60) 



which is well below the self-absorption frequency (see Medvedev & Loeb 2012 , for further details) 

~ (l0-VTC7enisMi?s/me)'/^^+'^i-2)/(^+^) 

~ 10^ - 10^ Hz, (61) 

where ax is the Thompson cross-section and we assumed a power-law energy distribution of elec- 
trons with index s with the nominal value of ~ 2.2. Note, however, that the peak frequency, i^g, is a 
very strong function of the NS surface field and the bubble size. The peak frequency also depends 
on time and, formally, exceeds 10^ Hz at At < 40 seconds before the merger. 

Although the spectral peak is self-absorbed, we can still use to calculate the non-absorbed 
part of spectrum, since P^, = Pu,max{'^ / '^s)^^'^^^''^'^ for the power-law distributed electrons. The 
spectral power at the peak (measured in erg s^-*^ Hz"-"^) is Pu,max{t) ^ P/vg = (c^meC^/Se)!?, 
where P is the total emitted power by a relativistic electron, P = (4/3)crT'C7g (iJ^/Svr). 

The unabsorbed observed spectral peak flux from a source located in our galaxy at the distance 
D = 10^2 cm (i.e., ~ 3 kpc) would be F^^^max = NePy^ma.^/i^'KD'^), where Ne = nejsMKhock is the 
total number of emitting electrons, hence 



F^,max(At) ~ 1^ (^y^shock^ISM j B 

~ 0.6 D^ie]i^fi3im-'/%sM,oRl%^t-'/' Jy- (62) 
The spectrum above the self- absorption frequency scales as 

F^{u, t) oc j,-(^-i)/2At-(3-5«)/8^ (63) 

For the nominal value of s = 2.2, we can estimate the observed flux at u = 10^ Hz as follows 

F,{At) ~ 0.01 Z?2"2'ee-'e°///lm-5nfsOM2 i^-^^z^g-O-eAt^- 1 Jy. (64) 
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Fig. 2. — Predicted spectral flux at 10^ Hz as a function of time. The NS-NS merger occurs at 



t = 0, hence At = -t in Eq. (64) 



This dependence is shown in Fig. [2j Note that this flux is very sensitive to the masses of the binary 
members, the system size (larger bubbles are fainter) and the strength of the surface fields of NS 
(hence binaries with magnetars are much brighter, though Rg for them is generally larger). 



7. Discussion and conclusions 



In this paper we elucidated one aspect of the question of what happens around a merging 
binary of compact objects — neutron stars, magnetars, black holes — with either one or both 
binary members being magnetized (i.e., we did not consider a double black hole binary). The binary 
dynamics is known to be determined by the gravitational wave emission, which happens within a 
finite time, hence the orbital period has, formally, a finite time singularity (neglecting general 
relativistic, finite size, and other effects). Rapid orbital motion of the objects' magnetospheres also 
lead to substantial electromagnetic (or Poynting flux) losses amounting to Eem ~ 10^^ ergs for 
NSs with the surface magnetic field of B ~ 10^^ gauss. This energy is extracted with relativistic, 
magnetized plasma, which forms an "electromagnetic bubble" of size Rg ~ 10^^ cm. The bubble is 
expanding at an accelerated rate, as the binary approaches the merger time, and drives a shock in 
the ambient medium. For the finite-time-singular law of the Poynting luminosity, we have obtained 
the analytical solutions for the bubble-|-shock evolution as a function of time at both early and late 
times. Using these solutions we were able to address observational signatures of such systems. We 
have found that they can be observed as faint radio sources. The spectral flux from a source about 
3 kpc away is estimated to be about a millijansky at ~ 10® Hz within At ~ 10^ seconds before the 
merger. The flux is increasing approximately as oc 1/At, so at At ~ 100 seconds before the merger, 
the flux is ~ 0.1 Jy. The observed flux is a strong function of the NS magnetic field oc 
and the shock radius oc Rj^- The expected radio signal for magnetars is detectable by existing 
radio observatories. Due to the large total electromagnetic energy release Eem ~ 10^'^ erg, the 
magnetar systems can be several tens times brighter than NS binaries. This estimate is, however. 
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rather uncertain because of the unknown shock/bubble size, which is expected to be about an order 
of magnitude larger than for NS binaries. 

To make observational predictions, we assumed the Poynting luminosity index to be p = 3/2. 
However, the actual value of p depends on the mechanism of the electromagnetic energy extraction 
and may differ from 3/2. This index can be determined from observations via simultaneous monitor- 
ing of the flux as a function of time and frequency ( Medvedev Sz Loeb||2012 ). If F,y{t) oc u^" (At)^* , 
then the energy injection index is 

p = (l-5/3, -2/3t)/(l-5/3,). (65) 
Interestingly, the predicted radio emission light-curve differs from other radio emission signatures 



of neutron stars and their binaries, e.g., the pulsar wind nebula emission (Piro & Kulkarni 2012) 



and the emission from the interaction of ejecta with the ISM (Nakar & Piran 2011). Particularly, 



our model predicts the very specific dependence of the flux on At, see Eq. (63), and the relation 



of the spectral and temporal indexes to the Poynting flux injection index, see Eq. (65). 



The scalings given in the paper are done for the shock co-moving time At, i.e., they do not 
include the finite light travel time from different parches of the spherical shock to the observer. 
When this effect is included, the radio signal from the shock of size ~ 10^^ cm will be spread 
over time ~ 2Rs/c, which is several hours. Thus we make a prediction that a short GRB should 
be accompanied by radio emission, which starts a few hours before the main event, thus being a 
precursor, and lasts for several hours after it. The question of what happens after the merger, 
goes beyond the scope of this paper. We can speculate that since GRBs form narrow jets, the 
nearly spherical bubble will not be affected substantially, except within the jet. After the merger, 
the Poynting luminosity vanishes, so the bubble pressure drops due to expansion and the forward 
shock eventually loses pressure support and starts to expand freely. Not too long after the merger, 
when/if the magnetic energy is still large, the magnetically-driven GRB scenario may be realized 
( [Lyutikov Sz Blandford 2003). Long after the merger, the shock should settle onto the Sedov-von 



Neumann- Taylor solution of a point-source explosion with energy E ~ 10^^ erg (neglecting other 
possible energy inputs), thus resembling a "micro supernova remnant". 

Several simplifying assumptions have been made in this study. In addition to the spherical 
symmetry, uniform external medium and neglect of general relativistic effects, there are a few others. 
Throughout the paper, we considered the shock to be strong, an assumption which may be wrong at 
early times (when the bubble expansion is slow because of low Poynting luminosity) and in relatively 
hot external medium, where the sound speed is comparable or larger then the bubble expansion 
speed — in the latter case, the shock will not form at all. We have also assumed that the shock is 
non-relativistic. This is a good assumption for NS binaries, but may be violated in magnetar binaries 
at late times. We have also assumed pressure equilibrium throughout the system, so the finiteness 
of the light travel time (which is roughly the time of pressure equilibration) has not been taken into 
account. Numerical simulations are needed to address such time-dependent evolution. One can 
expect, for example, the formation of a reverse shock propagating into the bubble. We have also 
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assumed that standard estimates of electromagnetic losses hold approximately true in the force- free 
magnetic configurations of orbiting companions. We based this assumption on numerical force-free 
simulations of pulsars, which generally confirm it. However, an accurate numerical modeling of both 
force-free structure and electron-positron cascades of orbiting magnetospheres is needed. Another 
interesting phenomenon here is that during the accelerated phase of the bubble-|-shock evolution, 
the bubble-shock interface can be prone to Rayleigh- Taylor instability, which can enhance mixing 
of the bubble material and the shocked ISM. Finally, we did not consider plasma dispersion cff'ects 
which may alter the travel time for the radio signal compared to the prompt and afterglow emission 
at other energies. 

We also note here on plasma processes at shock and details of particle acceleration, which 
impose additional constraints as follows. If the characteristic dynamical time of the system, which is 
^ At, is longer than the inverse collisional frequency f^^Jj then a collisional shock forms. Otherwise, 
when At < i/^jj, the shock is coUisionless. In this case, if At is (much) shorter then the Larmor 
frequency in the ambient field, the shock structure is not sensitive to the ISM field, but, instead, is 
determined by kinetic plasma instabilities (e.g., electrostatic Buneman or electromagnetic Weibel- 
like ones driven by particle anisotropy at the shock). The shortest associated timescale is the 
ion plasma time u^-^ lO^n^/^ s and, moreover, it takes about a hundred u^-^ seconds for an 
electrostatic shock [or uj^iv/c) seconds for a Weibel shock] to form and respond to changing 
conditions; it takes even longer for particle Fermi acceleration. Thus, the synchrotron shock model 
should be used with great caution for At as short as a fraction of a second or less. 

One of the authors (MVM) is grateful to Sriharsha Pothapragada for useful suggestions. This 

work was supported in part by DOE grant DE-FG02-07ER54940 and NSF grant AST-1209665 (for 
MVM), and NSF grant AST-0907890 and NASA grants NNX08AL43G and NNA09DB30A (for 
AL). 

REFERENCES 

Berger, E., 2011, New Astronomy Reviews, 55, 1 
Bildsten, L., k Cutler, C. 1992, ApJ, 400, 175 

Blinnikov, S.I., Novikov, I.D., Perevodchikova, T.V., Polnarev, A.G. 1984, Sov. Astron. Lett., 10, 
177 

Eichler, D., Livio, M., Piran, T., Schramm, D.N. 1989, Nature, 340, 126 

Etienne, Z.B., Liu, Y.L., Paschalidis, V., & Shapiro, S.L. 2012, Phys. Rev. D, 85, 064029 

Goldreich, P., & Julian, W.H., 1969, ApJ, 157, 869 

Gunn, J.E., & Ostriker, J.P., 1969, Nature, 221, 454 



- 18 - 



Lyutikov, M., & Blandford, R. 2003, preprint ( |astro-ph/0312347D 
McWilliams, S.T., & Levin, J. 2011, ApJ, 142:90 



Medvedev, M.V., & Loeb, A. 2012, preprint ( jarXiv:1212.0330D 
Nakar, E., 2007, Phys. Rep., 442,166 
Nakar, E., & Piran, T. 2011, Nature, 478, 82 
Piro, A.L. 2012, ApJ, 755, 80 



Piro, A.L., & Kulkarni, S.R. 2012, preprint ( |arXiv:1211.0547 ) 
Pshirkov, M.S. & Postnov, K.A. 2010, Astrophys. Space Sci., 330, 13 
Spitkovsky, S., 2006, ApJ, 648, L51 
Timokhin, A.N., 2006, MNRAS, 368, 1055 
Timokhin, A.N., 2010, MNRAS, 408, 2092 



This preprint was prepared witii tiie AAS lAT^rjX macros v5.2. 



